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We study rooted self avoiding polygons and self avoiding walks on deterministic fractal lattices 
of finite ramification index. Different sites on such lattices are not equivalent, and the number 
of rooted open walks W n (S), and rooted self-avoiding polygons P„(S) of n steps depend on the 
root S. We use exact recursion equations on the fractal to determine the generating functions for 
P„(S), and W„(S) for an arbitrary point S on the lattice. These are used to compute the averages 
< P n {S) >,< W n {S) >, < logP n (S') > and < log W n (S) > over different positions of S. We find 
that the connectivity constant /x, and the radius of gyration exponent v are the same for the annealed 
and quenched averages. However, < log P n (S) >~ nlog/i+ (a q — 2)logn, and < \ogW n {S) >~ 
nlog/i + (7 9 — l)logn, where the exponents a q and 7 9 take values different from the annealed 
case. These are expressed as the Lyapunov exponents of random product of finite-dimensional 
matrices. For the 3-simplex lattice, our numerical estimation gives a q ~ 0.72837 ± 0.00001; and 
7 9 ~ 1.37501 ± 0.00003, to be compared with the annealed values a a = 0.73421 and 7 a = 1.37522. 



I. INTRODUCTION 

Understanding the behavior of linear polymers in random media has been an important problem in statistical 
physics, both for reasons of theoretical interest, and applications. Calculation of quenched averages over the disorder 
is a very hard problem both analytically and computationally 1-2 . There have been many speculations and controversies 
regarding critical behavior of self-avoiding walks (SAWs) in random media (especially at the percolation threshold). 
Exact calculation of quenched averages has not been possible so far for any nontrivial case, and simulations are not 
easy and often give contradicting results for this problem. In this context, it seems useful to construct a toy model, 
where one can explicitly calculate the quenched and annealed averages, and see their difference. This is what we shall 
do in this paper, for the problem of linear polymers on a deterministic fractal. 

It is straight forward to calculate annealed averages of self-avoiding walks on deterministic fractals using real-space 
renormalization techniques^. One can explicitly write down a closed set of exact renormalization equations in a 
finite number of variables, so long as the fractal has a finite ramification index. Then, the eigenvalues of the linearized 
recursion equations near the fixed point of the renormalization transformation determine the critical exponents of the 
problem. In fact, a good deal of understanding of the complex behavior of polymers with additional interactions, e.g. 
self-interaction, or with a wall, or with other polymers, has been obtained by studying the corresponding analytically 
tractable problem on fractals 5 . 

It seems reasonable that the study of effect of inhomogeneities of the substrate would also be more tractable on 
fractal lattices. This is specially promising, as one does not need to introduce disorder in the problem from outside. 
The fractal lattices do not have translational symmetry, and hence a polymer living on a fractal lattice necessarily 
sees an inhomogeneous environment. Some regions of the lattice are better connected than others, and the local free 
density of the polymer per monomer in these regions would be lower than other parts. We want to understand the 
effect of presence of such regions on the large-scale structure and properties of the polymer. The main difference from 
the usual polymer in disordered medium problem to the case we study here is that the favorable and unfavorable 
. t-H ! regions are not randomly distributed over the lattice, but have a predetermined regular structure in the case of 
deterministic fractals. In this context, the annealed averages, which means averaging the partition function of the 
polymer over different positions, are appropriate in cases where the polymer can move freely over different parts of 
the lattice. The quenched average is averaging the logarithm of the partition function of the polymer over different 
positions, and would be appropriate where this freedom is not present. 

The annealed average for linear and branched polymers have been calculated exactly for many different fractal 
lattices 5 . But to our knowledge the quenched averages have not been calculated so far on any fractal lattice. In this 
paper, we use the recursive structure of fractals to calculate quenched averages for linear polymers on deterministic 
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fractals. We find that the connectivity constant \x, and the radius of gyration exponent v are the same for the annealed 
and quenched averages. The critical exponents for the quenched case can be expressed as the Lyapunov exponents 
for random product of finite-dimensional matrices. These can be estimated numerically efficiently by Monte Carlo 
methods, which we do for the illustrative case of 3-simplex fractal. 

The rest of the paper is organized as follows: In section 2 we define the 3-simplex lattice and specify the scheme we 
use to label the sites of the lattice. We also define the annealed and quenched averages precisely, and the generating 
functions for different quantities of interest. In section 3 we work out the recursion equations for generating functions 
of self-avoiding polygons(SAPs) and SAWs on 3-simplex. In section 4 we derive rigorous bounds on the number 
of rooted SAPs and SAWs, and prove that the connectivity constant /i, and the size exponent v are the same for 
the quenched and annealed averages. In section 5, we study the variation of the number of SAPs and SAWs with 
the position of the root on the 3-simplex lattice numerically. In section 6, we determine numerical values of the 
critical exponents in the quenched case by Monte Carlo determination of Lyapunov exponents for random products 
of matrices. 



II. PRELIMINARIES AND DEFINITIONS 

We will illustrate the general technique by working out explicitly the simple case of SAWs and SAPs on the 3-simplcx 
fractal. The treatment is easily generalized to other recursively defined fractals of finite ramification index. 

The 3-simplex graph is defined recursively as follows [Fig^: the graph of the first order triangle is a single vertex 
with 3 bonds. The (r + l)th order triangle is formed by joining graphs of three r-th order triangle by connecting a 
dangling bond of each to a dangling bond of the other rth order subgraphs. There is one dangling bond left in each 
graph and 3 bonds altogether. In general, the rth order graph will have 3 r_1 vertices and (3 r — 3)/2 internal bonds, 
and 3 boundary bonds. Clearly, the fractal dimension of this lattice is log 3/ log 2. 

We use a ternary base single integer to label different sites of an r-th order triangle. The labeling is explained 
in Fig. ^ A point on the r-th order triangle is labelled by a string of (r — 1) characters, e.g. 0122201 . . .. Each 
character takes one of three values 0,1 or 2. The leftmost character specifies in which of the three sub-triangles the 
point lies (0, 1 and 2 for the top, left and right sub-triangle respectively). The next character specifies placement 
in the (r — l)-th order sub-triangle, and so on. On an infinite lattice, specification of S requires an infinitely long 
string. In discussing the local neighborhood of a site, we only need to know the last few digits of S. We will denote 
by [S] r the substring consisting of the last r characters of the integer label of <S, and by s r the rth digit in the string 
S counted from the right. As an example, for the string S = 21. ...0112011, [5] 2 = 11, and S3 = 0. Also, in an obvious 
notation, [S] r = s r [S] r -i. 

The rooted polygons are polygons that pass through a given site (see Fig. ^1 called the root. Let P n (S) be the 
number of rooted polygons of perimeter n corresponding to root S. The generating function for rooted polygons for 
root S is defined as 



P(x;S) = ^2P n (S)x n 

n— 3 

Similarly, we define W n (S) as the number of open walks of length n whose one end-point is S, and the corresponding 
generating function W(x; S) by 



W(x;S) = J2Wn(S)x n (2) 

n=l 

For a given value of n, the values of P n (S) and W n {S) depend only on the last few digits in the ternary integer 
labeling S. For example, it is easy to check that if the last two digits of S are unequal, i.e. if s\ 7^ S2 we have 

P(x; S)=x 3 +x 6 + 3x 7 + 3x 8 + a; 9 + 0(x 12 ). (3) 

But if s\ = S2 

P(x;S) = x 3 + x 7 + 2x 8 + x 9 + 0{x 12 ) (4) 

In general, for any SAP configuration consisting of perimeter n, with n < 3 x 2 r , with r an integer, one can find 
an (r + l)-th order triangle graph, such that the polymer lies completely inside it. Then, the numbers P n (S) and 
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FIG. 1: The figure shows a 3 simplex graph of order 5, which is formed by joining 3 fourth order graphs. In general the rth 
order 3-simplex is formed by joining 3 graphs of (r — l)th order such that there is one dangling bond left in each subgraph. The 
thicker line shows a polygon of size 35 which passes through a particular site S. 



W n (S) depend only the relative position of S within this triangle, and hence only on [S] r . This then allows us to 
define averages of functions of P n (S) and W n (S). One assumes that S is equally likely to be any one of the 3 r sites 
in the (r + l)-th order triangle. So, for example, for n = 7, we have r = 2, and -FV(<S) depends only on [S]2- Out of 
9 possibilities for [S]2, 3 have s\ = 82, and 6 cases have s\ =/= S2- Thus, if S is chosen at random, using Eq. (3-4), 
we have Prob[P'j(S) = 3] = 2/3, and Prob[P-j(S) = 1] = 1/3. We shall use angular brackets () to denote averaging 
over different positions of the root S. This gives (P?) = 7/3, and (logfV) = flog 3. Other averages can be calculated 
similarly. 

We shall call the values (P n (S)}, and (W n (S)} as the annealed averages of P n (S) and W n (S), and define their 
generating functions 

00 

P( X ) = J2(Pn(S))x n (5) 
n— 3 



W(x) = ^(W n (S))x n (6) 

71=1 

We define the growth constant fi a , and the annealed exponents a a and -f a in terms of the behavior of (P n (S)) and 
(W n (S)) for large n: 

log(P n {S)) - nlog Ma + K-2)logn + 0(l) (7) 
log<W n (S)) = nlo gMa + ( 7a -l)logn + 0(l) (8) 

It was shown in ref.3 that fi a w 1.61803, a a w 0.73421 and -f a w 1.37522. 
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FIG. 2: Different restricted partition functions for the rth order triangle. Internal vertices inside the rth order triangle are not 
shown. 

For the quenched averages, the exponents a q and 7 q are defined by the condition that for large n 



We define the order of a polygon as the order of the smallest triangular subgraph that contains all the sites occupied 
by the polygon. For defining the size exponent v a and u q , it is sufficient to adopt the simple definition that the diameter 
of a polygon is 2 r if its order is r. One can then define the mean diameter of all polygons of perimeter n rooted at 
a given site S, as the average diameter, with all such polygons given equal weight. We define the quenched average 
mean diameter as the average over different positions of S of the average diameter of polygons of perimeter n rooted 
at S. The size-exponent v q for quenched averages is defined by the condition that the quenched average diameter 
varies as n Uq for large n. 

To define the annealed average for the diameter of SAP of perimeter n, we assign equal weight to all such loops 
within an sth order triangle with 2 s > n, and calculate the average diameter, ft is easy to see that the answer does 
not depend on s. Again, we define the annealed size exponent v a by the condition that the annealed average diameter 
varies as n" a for large n. 

The exponents for open walks can be defined similarly. We shall argue that the size exponents for open walks and 
polygons are the same, and further that \i q = /i a , and v q = v ai and hence simply write fj, and v without any subscript, 
if the distinction is unnecessary. 

Note that we have to first average over different positions of the root for a fixed n, and then let n tend to infinity 
to define v q . If we take the large n limit first, for a fixed position of the root, then even the convergence of large n 
limit of log[P n {S)fi~ n ]/ logn and log{W n (S)fi~ n ]/ logn is not obvious due to the irregular variation of log P„(S) and 
log W„(S) with n. The amount of fluctuations in different averages of observables over different positions of the root 
will be discussed in section 5. 



In this section, we briefly recapitulate the renormalization scheme for calculating the annealed averages used in rcf. 
3, and then adapt it for calculating properties of rooted walks. 

Consider one r-th order triangular subgraph of the infinite order graph. It is connected to the rest of the lattice by 
only three bonds. Our aim is to sum over different configurations of the SAW that lie within the subgraph, with a 
weight x for each step of the walk. These configurations can be divided into four classes, as shown in Fig. [21 and we 
define four restricted partition functions A^ r \ B^ r \ C^* and corresponding to these four classes. 

Here is the sum over all configurations of the walk within the r-th order triangle, that enters the triangle from 
a specified corner, and with one endpoint inside the triangle. B^ is the sum over all configurations of walk within 
the triangle that enters and leaves the triangle from specified corner vertices. CM and £>M are defined similarly (Fig. 
|2*|) . For any given value of r, flW^W^W^CW yfx are finite degree polynomials in x with non-negative integer 
coefficients. It is easy to see that the starting values of these variables are 



We note the generating function for P(x) is simply related to the generating function of unrooted polygons P n oroot(x) 
by the relation P(x) = x-^P noroot (x). The sum over all unrooted polygons of order (r + 1) within an (r + l)-th order 



<logPn(S)> 
(l0gW n (S)) 



nlogfiq + {a q - 2) logn + 0(1) 
n log M 9 + (7g-l) logn + 0(1) 



(9) 
(10) 



III. THE RENORMALIZATION EQUATIONS 



A™ = Vi, = x, C« = = 0. 



(11) 



B (r) 2 B (r)3 

FIG. 3: Two possible ways of getting a polymer configuration of type B 
triangle is , and number of sites in the (r + l)-th order triangle is 3 r , hence we get 

oo 

Pnoroot(x)=Y / ^ rB{r)3 - ( 12 ) 
r=l 

here we have suppressed the x dependence of B. In rest of the paper, we will suppress the x-dependence of variables 
A, B, C, D, to simplify notation, whenever the meaning is clear from the context. 
The sum over unrooted open walks can be expressed similarly 

oo 

Wnorootfi) = ^]3- r [3A (r)2 + 3B (r UW 2 + 3B (r)2 D (r) ], (13) 

r=l 

It is straight forward to write down the recursion equations for these weights A^ r+1 \ B^ r+1 \ C , ( r+1 ) and £)( r+1 ' in 
terms of the weights at order r. For example, Fig. [3] shows the only two possible ways one can construct a polymer 
configuration of type B. 



Similarly, we get 



\ C<r+V )-[ B M 2 ) \ C<r) ) (15j 

and 

D (r+i) _ A {r)2 + 2A {r)2 B [r) +AA {r) B {r) C {r) 

+65 (r) C (r)2 + D {r \2B^ + 3B (r)2 ) (16) 

We note that the equation i|14fl has one nontrivial fixed point B = B* = g" 1 ■ For starting value x < B* , the 
recursions give B^ -> as r tends to infinity, but for x > B* , B^ increases to infinity for large r. This gives 
the connectivity constant fi = 1/B* = (\/5 + l)/2, the golden mean. Linearizing the recursion equation about this 
nontrivial fixed point, we see that deviations from the fixed point value increases as £?( r+1 ) — B* w Ai(i?( r ) — B*), 
with Ai = 2 + ^- 2 . This then implies^ that v a = log2/logAi » 0.79862. 

The critical exponent 7 is determined in terms of the larger eigenvalue A2 of the 2x2 matrix in Eq. 15 evaluated 
at the non-trivial fixed point 3 , (A*,B*, C*) = (0, 1/fj,, 0). We get 



(17) 
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FIG. 4: Definition of weights for rooted graphs. In case of D, S can be any of the two end-points of the walk (shown by filled 
circles in the figure). 



The recursion equations for rooted SAP's and SAW's are constructed similarly. We define Bj (S) for s = 0, 1 and 
2 as the sum over walks on the r th order triangle that go through two corners of the triangle, visit the site S, avoiding 
the top, left and right corners for s — 0, 1 and 2 respectively (Fig^J. Since Ba T \s) depends on S only through [S] r _i, 
we shall write Bs(S) = Ba([S\ r —i). For any given r and S, these are finite degree polynomials in x. Similarly, we 
have to define three functions ([5] r _i) and Cj, ([5] r _i) for s = 0, 1, 2 instead of the single variables A, C and D 
for the unrooted problem, as the root breaks the symmetry between the corner sites of the r-th order triangle (Fig. 

mi. 

A site in the (r + l)-th order triangle is characterized by a string of r characters. Hence, a site characterized by 
string S at the r-th stage will be characterized by one of the strings OS*, IS* or 2S at (r + 1) th stage. 

We can now write the recursions similar to Eq. ^]for Bi r+1 \[S] r ), for s = 0, 1 and 2 in terms of weights of rth 
order graphs. For example, one can construct B^ 1 ^ (OS) in just one way as shown in Fig. [SJ The recursion relation 
for B { r+1) ([S] r ) will be 



4 r+1 )(0[S] r _ 1 ) = (sW) 2 5«([S] r _ 1 ) (18) 

Clearly, all the recursion relations for rooted polygons are linear in the rooted restricted partition functions of lower 
order Bs ([S\ r —i) , and can be written in the matrix form. With [S] r = s r [S] r —i, we have 



f B { ; +1 \[s] r ) 



B[ r+1 \[S} r ) 



\B { 2 r+1) ([S} r ) 



= M Sr I B[ r \[S\ 



(19) 



Bi r \[S], 



where 



( B 2 \ 

Mo = B b2 ; 

\ B 2 B J 



( B 0B 2 \ 
Mi = B 2 ; 
\ B 2 B J 
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M 2 



B B 2 
B 2 B 
OB 2 



(20) 



Here we have suppressed the superscript (r) of _B in the matrices Ado, M.\ and M. 2 . The generating function for 
rooted polygons, rooted at a site S is given by 



P(x,S) = Y j B£([S] r . 1 )B( 



r=l 



Similarly, we can write the recursion equations for A),' ' {[S] r -i) and G'i r \[S] r ~i) and D, 
to B^ {[S] r -i), with s = 0, 1, or 2. The A's and C's depend on each other and the recursion relations for them arc 



(21) 

^([S^r-i) defined analogous 



/ A {r+1] 

4 r+1 ' 

4 (r+i; 



([S]r) \ 
(lS]r) 
(lS]r) 
(iS]r) 
([S]r) 
([S]r) J 



where 



( ^([^r-l) \ 
At\[S]r-l) 
([5]r-l) 
CT([5]r-l) 



4 r - 

(r 



l) 



(22) 



£ = 



/ 1 B 2 B 2 \ 

B B 2 

B 2 B 

B 2 5 2 

B 2 

V0 B 2 ) 
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/BO B 2 
1 
B 2 B 

B 2 

\ 



\ 
B 2 B 2 

B 2 
B 2 
B 2 J 



C 2 



( B 


_B 2 














B 2 


B 
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B 2 


B 2 














B 2 




















B 2 





\o 





B 2 








B 2 



(23) 



Here again we have suppressed the superscript (r) on the _B's. One can now write similar recursion for D^ r \[S] r ) 
also. [Like the annealed case£, these variables are not needed for the determination of critical exponents, though they 
are needed to determine W rl (5').] Here we write down the recursions for D^ +1 \[S] r ). similar relations will hold for 



D{ r+1 \[S] r 



and D 



(r+l) 



(lS]r). 



^ r+1) (0[S] r _i) = {AB + 2BC){A l +A 2 ) + B 2 D 

^ r+1) (l[5] r _i) = (A + AB)A 1 + (AB + 2BC)(C 1 +C 2 ) 
+ D a B + D 2 B 2 + 2BCC 

{A + AB)A 2 + (AB + 2BC)(d + C 2 ) 



D, 



(r+l) 



(2[S]r-l) 



-D B + D\B + 2BCC 



(24) 



where we have suppressed the ([S] r _i) dependence, and the superscripts (r) in the terms on the right hand side. 

We can also express the open walks generating function W(x, S) in terms of these restricted partition functions. It 
is easy to show that 



W(x;S) = J2 



A a > ([S] r -i)(A + AB) + D Sr ([S] r -i)B 2 

S ^S r 



+ AA Sm ([S] m -i])+0(x 2m ") 
In the limit of large to, the last term can be dropped, and we get 



(25) 



W(x;S)= Y 



K' ({S] r -i)(A + AB) + D Sr ([S] r ^)B 2 



s ^S r 



(26) 



For order r, there are 3 r_1 different choices for [S] r _i, and 3 choices of s each of the A S ,B S , C s and D s , a total of 
4.3 r polynomials. The generating function for walks rooted at one end W(x, S) can be easily written in terms of these 
generating functions. To determine P(x, S) and W(x, S) to given order n in x, we need to determine these only upto 
a finite order r. Thus, we have an efficient algorithm to explicitly determine these polynomials to any given order, 
and also to calculate the different averages over different sites. The explicit enumeration of rooted walks on fractals 
was studied by Reisi and Ordemann et. al£, but the sizes they could reach were much smaller. We have obtained 
exact values of SAWs upto length 128 and SAPs upto 768 for all sites of the fractal, using our recursion equations 
(see Table. 1) 
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IV. DETERMINATION OF THE GROWTH CONSTANT 

We note that the recursion equation of (Eq. 14) does not involve the rooted variables B s r (S), but not vice 
versa. If we start with a value x < B* = 1/fJ., then for large r, flW tends to zero. This would make the rooted 

variables Bi r \S) also tend to zero for large r. However, if x > B* , then B^ will diverge for large r, and this would 

make Bs r \s) also diverge for all S. Thus we conclude that for all S, P{x, S) and W(x, S) converge if x < and 
diverge if x > l//i. Hence it follows that the critical growth constant is /i, for all S. 

Consider now polygons or walks of finite length n, with n large. We think of x as the fugacity variable for each 
step. Then large n corresponds to initial value x very near l/fi, say x = 1/fi — S where 5 scales as 1/n. Under 
renormalization, the value of B* — B^ increases, and there is a value tq such that for r — r$ 1, flW ps 0, and for 

ro — r l,£>( r ) ps B* . Then, the recursion equations 1191 and 1221 imply that the rooted partition functions Bs T \S) 
and Ci r \S) also become very small for r > r , while As tends to a finite value as r tends to infinity. The value of 
r increases as 6 is decreased as r ~ ifg^p 1 an d the diameter of polymer increases as 2 r ° ~ 8~ Va . As 8 ~ 1/n, this 
implies that the average size of polymer increases as n v ', independent of S, and hence u q = v a . 

An alternative proof of the assertion that the growth constant [i is the same for all sites is provided by setting up 
upper and lower bounds for P n (S) and W n (S), which have the same exponential growth. 

We first obtain an upper bound on P n (S). If we ignore the constraint that the walk has to pass through S, we 
get an upper bound on the number of such walks. For example, for walks contributing to B^ (£) we can write 

Ba(x,S) < B^(x) where a = 0,1,2 and the inequality between polynomials is understood to imply inequality for 
the coefficient of each power of x. This implies that for all sites S, 

00 

P(x,S)<Y / B {r)3 (27) 

r=l 

If we write the function giving the upper bound in right-hand-side as U(x), then U(x) satisfies the equation 

U(x) = x 3 + U{x 2 +x 3 ) (28) 

This functional equation again has the fixed point at x* — 1/fJ-, and linear analysis near the fixed point shows that 
U(x) diverges as — log(l — xfi) as x tends to l//z from below. This implies that the coefficient of x n in the Taylor 
expansion of U(x) varies as \x n jn for large n. Thus, 

P n {S) < K^ n /n, for all n, and all S, (29) 

where K is some constant. 

We now obtain a lower bound for P n (S). Let 5*0 be the site whose label is a string with all digits 0. This is the the 
topmost site of the triangular graph. For such a site Bq T \x; Sq) satisfies the following recursion 



Bi r \s ) = (B( r -VfBt 1) (S ) (30) 

Clearly, for all sites 5, we have B { 1} (S) > B^\s ). Then, using Ea. (|18|) and (|19|l . by mathematical induction, we 
see that for all r, s and S 

Bi r \S)>B^(S Q ), for S = 0,l,2. (31) 

This implies for any site S 

P(x,S) > P(x,S )=J2 B{r)2B o\s ) (32) 
r=l 

Clearly, the lower bound is actually attained for S — So. If we take P(x,Sa) — xL(x), then L(x) satisfies the 
following equation 



L(x) = x 2 + x 2 L(x 2 + x 3 ) 



(33) 
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Assuming that L(x) near x = l/fj, has a singular expansion of the form L(l//x — 5) = L(l//x) — ^STir, where K is 
some constant, we get 

6 = 21og^/log(2 + /i- 2 ) = 0.92717 (34) 

and hence 

Pn(S) > Ku^n-'- 1 (35) 
where K\ is a constant. Hence for any S, we have proved 

ifiM™™- 6 - 1 < P n (S) < Kfi n /n (36) 
Thus, in the limit of large n, for all sites S, we must have 

n m i2i^) =/i (37) 

We also get the nontrivial bounds 

1 - b < a q < 1 (38) 

Similarly it is easy to prove upper and lower bounds for open walks. A simple lower bound is provided by the 
inequality W n (S) > P n (S), for all S. An upper bound to (log W n (S)) is provided by \og(W n (S)} . But the latter is 
known to vary as nlog/i + — 1) logn for large n. Hence we get 

lim ± (log W n (S)) =/i. (39) 

n—>oo ft 

V. VARIATION OF P n (S) AND W n (S) WITH S AND n 

We use the recursion equations (19-25) to calculate the values of P n (S) and W n (S) for different choices of the root 
S. For an SAP of order r, the minimum perimeter is 3 x 2 r ~ 2 . We could study these recursions upto r = 9 and 
hence took into account all self avoiding polygons upto sizes n — 3 x 2 8 = 768. A 9-th order triangle has 3 8 = 6561 
vertices, we calculated the different polynomials A S (S), B S (S),C S (S) and D S (S) for s = 0, 1, 2 and 3 8 possible values 
of S, keeping all terms up to order x 768 in each polynomial. The coefficients for large order become rather large. For 
example, the number of polygons of size 768 for the topmost site of the fractal graph is « 2.5 x 10 154 . We used the 
symbol manipulation software Mathematica 6 , which allows one to work with integers of arbitrary size. These then 
are used to calculate various averages. 

In Fig. El we have shown the variation of P n (S)n~ n with n for some selected values of S. We see clearly roughly 
log-periodic variation in these numbers. We have also plotted the upper and lower bounds on P n {S) derived [Eq. l|27l 
and OH)]. As was argued there, there is a site So which saturates the lower bound. However, there is no single site 
that saturates the upper bound. Most of the polygons having a given value of perimeter n, have a particular order r, 
and have to pass through the three bonds joining the (r — l)-order subgraphs. The six sites that are at the ends of 
these bonds clearly maximize P n (S). The sites change if n changes to correspond to polygons with one higher order. 

In Fig. we show the variation of exactly calculated values of log < P n (S) > and < logP n (S) > with n. We also 
show the maximum and minimum values of logP n {S) attained, as a function of n, as S takes all possible values. Note 
that though we could study sizes upto 768 exactly, it is difficult to estimate a a and a q even to two digit precision 
from the data using standard series extrapolation techniques because of log-periodic oscillations. 

The log-periodic variation of P n is easy to understand qualitatively^*^. All SAP's of order r in a given (r + l)th 
order triangle have to pass through the three constriction points where the constituting r-th order triangles are joined. 
There is a natural length n ~ C\\, where C is some constant, for a polymer that does this. If the length is somewhat 
smaller than this value, the polymer is a bit stretched, and has a lower entropy. If it is a higher by a factor 1.5 or 
so, it loses entropy as many monomers have to squeeze in the same space. If the value of n increases by a factor Ai, 
then same thing happens at a higher order triangle. The sites which maximize P n (S) for a given n are the at the 
constriction points of the the corresponding r-th order triangles. As logn changes, these points also change. The 
fractional number of points S for which P n (S) attains its maximum value clearly varies as 3 _r , where r is the order 
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FIG. 6: Figure shows H n (S) = logP n (S)fj,~ n plotted as a function of n for five different values of S. These S are chosen near 
the constriction points and hence they maximize P n (S) for a given n. Here Si = 00000211, S2 = 00002111, £3 = 00021111,54 = 
00211111 and S5 — 02111111. Each of these values of S attains the maximum for some range of n. The points, defined only 
for integer n, have been joined by straight-line segments as an aid to the eye. 

of a typical loop of perimeter n. Using n ~ R we see that the density of points varies as R-° ~ 71-°/", where 
d = log3/log2 is the fractal dimension of the lattice. 

The log-periodic oscillations for open walks are much smaller in magnitude than in case of polygons. Here also, 
the number of open walks or order r is greater if the root is near the corner point of an r th order triangle, as then it 
gets more space to explore and hence more entropy. We find the amplitude of oscillations in H n is approximately 100 
times smaller than for the corresponding quantity for closed polygons shown in Fig. and [71 In Fig. [S]we show the 
exact values quenched and annealed averages for walks upto size 128. The difference in the two averages in the case 
of walks is smaller than for the polygons. 

In Figs. and we have shown the density plots for SAPs on fractal lattices for n — 400 and n = 600, The 
different sites are shown by different colours depending on the value of logP n at that site. The blue regions represent 
sites with largest values of P n (S). We clearly see that depending on value of n some sites are more favoured than 
others. We note that the difference between maximum and minimum values of log P n (S) is more than 10 for n = 400, 
but only about half of this value for n = 600. 

VI. CALCULATION OF EXPONENTS FOR QUENCHED AVERAGES 

We note that to a very good approximation, all loops of a given perimeter n have the same order , say ro, where 7*0 
is an integer approximately equal to ■ The contribution of these loops to P(x; S) then comes mostly from the 

single term T ro — Bs^' ) ([5 , ] ro -i)i? < - ro ' 2 corresponding to r = ro in the equation (|21(l . Also, if x = 1/fi = B* , then the 
contribution is P n (S)[i~ n , which is a slowly varying function of n. The contribution of nearby values of n to the ro-th 
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FIG. 7: A plot of average value of logP n (S)n~ n (H n ) as a function of n. The uppermost and the lowermost curves are 
theoretically derived upper and lower bounds to this number over different positions. The dashed and dotted lines show the 
annealed and quenched average value respectively. 

term will be comparable. The number of such terms that contribute to the ro-th term is of order n. Thus we have, 

P n (S)fi- n ~ ~T ro (S) (40) 
For x = B*, the matrices Mo, Mi, M2 become independent of r, and T ro is of the form 

T ro = (v 1 \M. r _ 1 ...M at M S3 M S2 M si \v 2 ) (41) 

Here {v\\ and (1)2) are specific 3-dimensional bra- and ket vectors. Then from the general theory of random product 
of matricesii*i2iiii ; it follows that the probability distribution of T ra tends to a log- normal distribution for large ro, 
and 

lim — logT ro =logAi (42) 

r ^oo ro 

where log Ai is the largest Lyapunov exponent for the random product of matrices. This also implies that the variance 
of logP„(S') will be proportional to logn. This is much less than some positive power of n, which is the expected 
behaviour for the usual polymer in the random medium problem. For example, for the directed polymer in a random 
medium, the variance of free energy of an n monomer chain varies as n 2 / 3 . This is due to the fact that in the 
deterministic fractal case studied here, the favourable and unfavourable regions are very evenly distributed. 
It then follows that (log[P n (5) M - ™]) tends to r logAi — logn. Putting r ~ (logn)/(logAi), we get 



a q = 1 +logAi/logAi 



(43) 




It is straightforward to estimate Aj numerically. We start with an arbitrary initial 3-dimensional vector \vq) of 
positive elements with sum 1, and evolve it randomly by the rule 

|t> J+ i) = a j M j \v j ) (44) 

where Mj is randomly chosen to be one of three matrices Mo, Mi, M2 and aj is a multiplying factor chosen so that 
the sum of elements of new vector is again 1 . The Mj are 



/ B* 2 \ 
M = B* B* 2 ; 
\0 B* 2 B* J 



Mi = 



B* 





J5* 





B* 2 





B* 2 








Ma = 





B* 2 


















B* 2 



(45) 



where = (V5-l)/2. 

We iterate ea!44l many times, and estimate logAi by logAi 
°"/V jmax, where cr is the rms fluctuation of logaj. Values of j„ 
jmax = 10 9 , we find a « 0.45 and 



^j™ 1 logaj. The error in logAi decreases as 
require less than a minute of CPU. For 



10 8 



logAj 



-0.23575 ±0.00001 



(46) 




FIG. 9: The plot shows the log(P n (S)) for n = 400 for all sites for fractal lattice upto r = 9 generations. 




FIG. 10: The plot shows the log(P n (S)) for n — 600 for all sites for fractal lattice upto r = 9 generations. 



16 



Note that this corresponds to a lattice with 3 Jmaa ' sites. This gives a q — 0.72837 ± 0.00001. 

A similar calculation for the exponent j q can also be done. To calculate the exponent j q , we note that again the 
most contribution to W(x; S) will come from terms with r as ro. At x = B* , C$,C\ and £2 become independent of r. 
If 



C/ ro = (u 1 \C Sr _ 1 ...C S4 C S3 C S2 C Sl \u 2 ) (47) 

where < U\\ and \u 2 > are 6-dimensional bra- and ket vectors and if log A 2 is the largest Lypanouv exponent for the 
random product of matrices, 



lim -log[/ ro =logA 2 (48) 

ro— »oo ro 

Hence, log^ r) (S*) and logCf \S) will increase as rlogA2, for almost all S, for r < ro. For r > ro, A^p w A^"" 1 
and C> « 0. Therefore it follows that the leading order contribution to equation (29) will come from r — tq term. 
Hence < logW n (S)^~ n > tends to /o<7(A2A 2 )r — log{n) and we get 



7 9 = (49) 
log(Ai) 

We find for j max = 10 9 , a s» 1 and A 2 = 1.04845 ± 0.00003. Substituting in equation above we get 7 g = 
1.37501 ±0.00003. 
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n 


log < P n > —nlogfj, 


< logPn > —nlogfj.) 


log < W n > —nlogfj, 


< logW„ > —nlogfj, 


3 


-1.4431783323 


-1.4442955737 


0.3481239636 


0.3481239636 


7 


-2.5207277722 


-2.6372273881 


0.6327712915 


0.5825518549 


8 


-2.8684082046 


-2.8875491290 


0.6648211368 


0.6276714199 


9 


-4.3304492827 


-4.3328867211 


0.6797287773 


0.6707214751 


15 


-4.0919348312 


-4.5179887766 


0.9155729406 


0.8693414702 


16 


-3.8277310472 


-4.0586132570 


0.9370278020 


0.8829866367 


17 


-3.7374926267 


-3.8857022562 


0.9544167885 


0.8994169039 


18 


-3.8145473883 


-3.9148386256 


0.9693553123 


0.9204691093 


19 


-4.0133264968 


-4.0828080387 


0.9815329752 


0.9429045509 


20 


-4.3366352924 


-4.3875115144 


0.9920263470 


0.9648015516 


40 


-4.8957208753 


-5.0810665168 


1.2742813939 


1.2185712350 


80 


-6.6481626291 


-7.2828262546 


1.5384807521 


1.4965133405 


120 


-7.0076576021 


-7.0786974211 


1.6575701703 


1.6393667050 


160 


-9.5287171328 


-10.9169115887 






320 


-9.6664681509 


-9.7339866481 






640 


-8.6892637799 


-8.7886721125 






768 


-10.8851067629 


-10.9571396026 







TABLE I: Values of quenched and annealed averages for SAP and SAW for some representative values of 



